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Abstract 

Background:19PM ion channels are involved in nociception and are also excited by pungent odorous substances. Based 
on reported associations of TRPAl genetics with increased sensitivity to thermal pain stimuli, we therefore hypothesized 
that this association also exists for increased olfactory sensitivity. 

A7ef/jO£/5; Olfactory function and nociception was compared between carriers (n = 38) and non-carriers {n = 43) of TRPAl 
variant rsl 1988795 G>A, a variant known to enhance cold pain perception. Olfactory function was quantified by assessing 
the odor threshold, odor discrimination and odor identification, and by applying 200-ms pulses of H2S intranasal. 
Nociception was assessed by measuring pain thresholds to experimental nociceptive stimuli (blunt pressure, electrical 
stimuli, cold and heat stimuli, and 200-ms intranasal pulses of CO2). 

/fe5i//f5; Among the 11 subjects with moderate hyposmia, carriers of the minor A allele {n = 2) were underrepresented (34 
carriers among the 70 normosmic subjects; p = 0.049). Moreover, carriers of the A allele discriminated odors significantly 
better than non-carriers (13.1 ±1.5 versus 12.3±1.6 correct discriminations) and indicated a higher intensity of the H2S 
stimuli (29.2±13.2 versus 21 ±12.8 mm VAS, p = 0.006), which, however, could not be excluded to have involved a 
trigeminal component during stimulation. Finally, the increased sensitivity to thermal pain could be reproduced. 

Condusions:J\r\e findings are in line with a previous association of a human TRPAl variant with nociceptive parameters and 
extend the association to the perception of odorants. However, this addresses mainly those stimulants that involve a 
trigeminal component whereas a pure olfactory effect may remain disputable. Nevertheless, findings suggest that future 
TRPAl modulating drugs may modify the perception of odorants. 
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Introduction 

Genetic mutations in ion channels encompass a diverse range of 
pathological conditions including extreme pain conditions such as 
total insensitivity to pain or exaggerated paroxysmal pain [!]■ A 
rare gain-of-function mutation in TRPAl transient receptor 
potential cation channels causes the familial episodic pain 
syndrome [2]. TRPAl is expressed on nociceptive neurons [3,4]. 
They are gated by noxious cold (<17°C [4]), although this is 
species-specific [5] and can be switched to heat sensitivity by 
changes in N-terminal ankyrin repeat-rich domain [6]. 

Due to their role in pain perception, TRPAl channels are 
reckoned as a promising novel target for new analgesics [7]. 
However, the function of TRPAl may exceed the somatosensory 
system. In fact TRPAl channels are expressed in the olfactory 
epithelium [8] and the olfactory bulb [9] of mice. This raises the 



possibility of a specific involvement of TRPAl in olfaction. In 
addition, among activators of TRPAl channels are also pungent 
odorous chemicals [7,10]. This further suggests a role of TRPAl 
in perception of odors with a trigeminal component. This is also 
supported by in vitro observations of TRPA 1 activation by odors 
[11]- 

To address an involvement of TRPAl in human olfaction, the 
present study employed a genetic approach based on the 
hypothesis that positive associations should then be multisensory. 
Specifically, gain-of-function TRPAl variants that increase the 
sensitivity to pain stimuli should also increase the sensitivity to 
odorous stimuli. A suitable variant accommodating this approach 
was found with the TRPAl rsl 1988795 G>A single nucleotide 
polymorphism. It has been associated previously with an increased 
sensitivity to nociceptive stimuli [12], however, without molecular 
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proof yet. The present work pursued the hypothesis that this 
enhanced sensitivity will extend to olfactory stimuli. 

Methods 

Subjects and Design 

The study followed the Declaration of Helsinki and was 
approved by the Ethics Committee of the Goethe-University 
Frankfurt am Main, Germany. Informed written consent from 
each participating subject had been obtained. The study largely 
paralleled a similar assc-ssment of a linkage between olfaction and 
nociception via a Na^l.7 sodium channel genotype [13j, without 
overlapping subjects. It was double blind with respect to the 
subjects' TRPAl genotype. The subjects' actual health was assessed 
by medical history, physical examination including vital signs, and 
routine clinical laboratory' test results. Exclusion criteria were a 
current clinical condition affecting pain or olfaction, any other 
actual diseases and drug intake within a week except oral 
anticonceptionals. 

Power calculations based on published data [14,15] had resulted 
in a sample size of 80. Due to the enrolment of three non- 
Caucasian subjects and one additional Caucasian, 84 subjects were 
enrolled, which was covered by the Ethics approval. The cohort 
submitted to analysis consisted of a random sample of 38 unrelated 
healthy men and 43 unrelated healthy women (age 20-30 years, 
mean ± standard deviation 24. 1 ± 2 years, 1 3 smokers) of 
Caucasian ethnicity by self-assignment. The minor rs 1198879,5 
A allele was carried by 36 subjects (allelic frcrjucncy 26.5%), was 
equally distributed among genders (cross-tabulation: likelihood 
ratio = 0.344, p = 0.38) and the distribution of 45 non-carriers, 29 
heterozygous and 7 homozygous carriers corresponded to the 
Hardy-Weinberg equilibrium (x^ goodness of fit test: p = 0.54). 
The tests were done in the succession pressure pain, heat pain, 
cold pain, application of menthol, olfactory testing (Sniffn' sticks), 
cold menthol pain, electrical pain, and intensity ratings of CO2 
and H2S intranasal stimuU (for method details, see below). 

Genotyping 

Genomic DNA was extracted from 200 ^,1 blood using the EZl 
DNA Blood 200 nl Kit on a BioRobot EZl Workstation (Qiagen, 
Hilden, Germany). Genot>ping for TRPAl rs 11988795 G>A was 
done after the end of the study from genomic DNA by means of 
Pyrosequencing assays on a PSQ,96 MA System (Qiagen, Hilden, 
Germany) using the PyroMark Gold Q96 Reagents set (Qiagen, 
Hilden, Germany). PGR reactions were performed in a 25 |J,1 assay 
volume on a Mastercycler ep gradient S instrument (Eppendorf, 
Hamburg, Germany), using the HotStar plus Taq Polymerase 
system (Qiagen, Hilden, Germany) and SNP-specific PGR primers 
(forward primer: 5'-TAAGTGAGCCAAGTTCAGATCAGA-3' 
and reverse primer: 5'-biotin-TTTCACAGAAAGT- 
GAGGTGTTGTA-3'). The PGR was done with an initial 
denaturation step for 5 min at 95°C, 50 cycles with a 30 second 
denaturation step at 95°C, an annealing step at 45°C for 30 
seconds and an elongation step at 72°C for 30 seconds, followed by 
a final elongation step at 72°C for 5 min. The PGR product 
(25 (J.1) was used in the Pyrosequencing analysis as described 
previously [16] with the sequencing primer 5'- 
TGATCCTTCTTTTCTCAGTA-3'. Three samples of each 
genotype were also sequenced by an external provider (EGG 
GmbH, Berlin, Germany), using the conventional capillary 
sequencing method [17] on a ABI 3730 XL device (AppUed 
Biosystems, Darmstadt, Germany), and implemented as positive 
controls during Pyrosequencing. 



Assessment of Olfactory and Nociceptive Function 

Olfactory testing. The olfactory test was based on felt-tip 
pens that contained a solution of an odorant instead of lirjuid dye 
("Sniffin' Sticks": Burghart, Wedcl, Germany [18]). The pen's cap 
was removed by the experimenter for approximately 3 s and the 
pen's tip was placed 1-2 cm in front of the nostrils, in the case of 
triplet pen presentation at an interval of approximately 3 s. Three 
main components of olfactory function were assessed birhinally, 
namely the perception of odors at low concentrations, which is the 
odor threshold, the ability to recall an odor and name it, which is 
odor identification, and the distinction of different smells, which is 
the ability of odor discrimination. 

Odor thresholds. Were obtained for the rose-like odor 
phenylethylethanol, which was presented in 16 successive 1:2 
dilution steps starting from a 4% solution (i.e., higher dilution steps 
mean less concentrated PEA). Using a three-alternative forced- 
choice task (3-AFC) and a staircase paradigm starting at low 
phenylethylethanol concentrations, one pen with the odorant and 
two blanks were presented at each dilution step. Two successive 
correct or one incorrect identification triggered the reversal of the 
staircase. The odor threshold was the mean of the last four out of 
seven staircase reversals. 

Odor discrimination. Was determined with 16 triplets of 
pens, two containing the same odorant and the third a different, 
"target" one (i.e., (target/non-target) butanol/2-phenyl ethanol, 
isoamylacetate/anethole, anethole/eugenol, limonene/fenchone, 
(— )carvone/(+)carv()ne, eugenol/cinnamon aldehyde, dihydroro- 
senoxide/ menthol, acetaldehyde/isoamylacetate, citronellal/linal- 
ool, pridine/Umonene, limonene/ citronellal, eucalyptol/ dipyridyl, 
dipyridyl/ cyclopentadecanoate, butanol/fenchone, octylacetate/ 
cinnamon aldehyde, carvone/acetaldehyde). Using a 3-AFG 
paradigm, the subject's task was to identify the target stick. 

Odor identification. Was determined with 16 odors (i.e., 
orange, leather, cinnamon, peppermint, banana, lemon, liquorice, 
turpentine, garlic, coffee, apple, clove, pineapple, rose, anise and 
fish) using a 4-AFC task with presentation of a list of four 
descriptors for each pen. From the three olfactory subtests, 
normosmia, i.e., normal olfactory function, was established by 
calculating a composite "TDI score" ("Threshold Discrimination 
Identification") as the sum of the scores from the three subtests 
[19]. Normosmia is observed at TDI scores >30.5 whereas 
pathologic olfactory function is indicated by TDI S30.5, with the 
separation of hyposmia (30.5STDI>15.5) from functional anos- 
mia at TDI <15.5 [14]. 

In addition to the olfactory test, suprathreshold olfactory stimuli 
were delivered to the nasal mucosa by means of an olfactometer 
(OM/2, Burghart Messtechnik GmbH, Wedel, Germany). Spe- 
cifically, a Teflon tube (outer diameter 0.5 cm) was introduced 
approximately 1 cm into the right nostril. Through this tube, 
pulses of H2S (n = 1 8 stimuli, length 200 ms, concentration 
5 ppm) were delivered directly to the mucosa [20] at 36-44 s 
intervals by embedding them in a constantly flowing airstream (8 
1/min) at controlled temperature (36.5°C) and humidity (80% 
relative humidity). During this experiment the subjects were 
comfortably seated in an air-conditioned and visually shielded 
room; acoustic shielding was achieved with white noise (50 dB 
SPE) delivered via headphones and subjects observed a special 
breathing technique that avoids airflow in the nasal cavity 
(velopharyngeal closure [21]) to ensure correct stimulus delivery. 
The subjects performed a tracking task on a computer screen [22] 
where within 2.5 s after each stimulus a visual analog scale (VAS) 
was displayed to query ratings of the stimulus' smell intensity ("no 
smell" to "maximum smell"). The median of the estimates of 18 
single stimuli was submitted to statistical analysis. 
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Experimental pain testing. Pain thresholds to blunt 
pressure were obtained using a pressure algometer with a 
circular and flat probe of 1 cm diameter (Commander Algometer, 
JTECH Medical, Midvale, Utah). It was placed perpendicularly 
onto the mid-phalanx of the right middle finger. The pressure was 
increased at a rate of approximately 9 N/ cm^ per second until the 
subject indicated pain. The increase in pressure was controlled 
manually by the investigator and stopped once the subject 
reported pain. The procedure was repeated five times at intervals 
of 30 s and the pain threshold was the median of the five 
measurements. 

Pain thresholds to electrical stimuli were obtained using a 
constant current device (Ncurometer® CPT, Neurotron Inc., 
Baltimore, MD) that delivered sine-wave stimuli at 5 Hz apphed 
via two gold electrodes placed on the medial and lateral side of the 
mid-phalanx (middle finger of the right hand as default-testing 
site). Their intensity was increased from 0 to 20 mA by 0.2 niA/ s. 
During the test, subjects kept a button continuously pressed until 
they felt pain and interrupted the current by releasing the button. 
Measurements were repeated five times at intervals of 30 s and the 
median of the electrical currents at which subjects had released the 
button was the pain threshold. 

Pain thresholds to cold stimuli were obtained using a 
Thermal Sensory Analyzer (Medoc Advanced Medical Systems 
Ltd., Ramat Yishai, Israel). Cold stimuh were applied with the 
3x3 cm^ thermode placed on a skin area of the left or right 
(randomized) volar forearm. The temperature was lowered from 
32°C to 0°C by l°C/s. The subject pressed a button at the first 
sensation of pain. Measurements were repeated five times and the 
median of these measurements was the threshold. Subsequently, a 
plaster soaked with menthol (2 ml of a 40% menthol cUssolved in 
90% ethanol) was applied onto the skin area for 30 min and cold 
pain measurements were repeated (cold/ menthol pain thresholds). 

Pain thresholds to heat stimuli were obtained using the same 
device as for cold stimulation, placing the 3x3 cm thermode on 
the contralateral volar forearm (to forearm used for cold pain 
measurements). Temperature was continuously increased from 
32-52. 5°C by 0.3°C/s and the subject pressed a button at the first 
sensation of pain. As previous experiments had shown that three 
heat applications are needed to reach a plateau in the thresholds, 
heat applications were repeated eight times and the median of last 
five responses was used as heat pain threshold. 

Finally, the perceived intensity of chemical pain was queried by 
using CO2 (60% v/v, n=18 stimuli) instead of HjS as 
olfactometer stimulus applied to the right nostril. CO^ is converted 
into bicarbonate and protons by the enzyme carboaiihydrase [23], 
it evokes a short stinging pain sensation [24] due to excitation of 
trigeminal nociceptors [25,26]. In analogy to the HjS-stimuli, the 
CO2 pulses were rated for their painfulness using visual analogue 
scales ("no pain" to "maximum pain"). The median of the 
estimates of 18 single stimuli was submitted to statistical analysis. 

Statistics 

Olfactory test results (threshold discrimination, identification 
and intensity estimates of H2S stimuli) were standardized by z- 
transformation and submitted to single analysis of variance for 
repeated measures (rm-ANO\'A) using "odor task" as within- 
subject factor (4 levels) and "genotype" as between-subjects factor. 
In addition, olfactory test results were compared separately 
between genotype groups by means of t-tests, without z- 
transformation. To accommodate the complexity of pain pheno- 
types by clustering as proposed previously [27-29], the otherwise 
similar analyses of the pain data included the additional factor 
"pain cluster", which was obtained by Ward hierarchical 



clustering with the squared Euclidian distance and subsequent 
description by classification and regression tree (CART) classifiers. 
Specifically, rm-ANOVA was calculated with z-transformed pain 

parameters (n = 6, cold and menthol pain thresholds multiplied 
with — 1 to obtain similarly directed measures), using "pain 
measure" as within-subject factor (6 levels) and "genotype", and 
"pain cluster" as between-subjects factors. To check for gender 
diflerences in the genotypic effects, "gender" was introduced in a 
further combined analysis, however, at an exploratory level as the 
study was not powered for gender differences. Moreover, while the 
original study reporting a function of the TRPAl rsl 1988795 G> 
A variant also reported interactions in the modulation of pain 
between gender and several variants in other genes [30], gender 
interaction was not reported for TRPAl variants and therefore not 
a major factor to be introduced in the present study design. 
Further statistical analyses included cross-tabulations, rank corre- 
lations (Spearman's p) and effect size calculations using Cohen's d 
(values of 0.2, 0.5 and 0.8 being indicative of a small, medium and 
large effect size, respectively [31]). A gene dose effect was not 
assessed for reasons of statistical power (nA,\ < ii,,icr., = 7). Statistics 
were done using SPSS (version 21 for Linux, SPSS Inc., Chic;ago, 
USA) or R (version 2.14.1 for Linux; http://CRAN.R-project. 
org/); the a-level was set at 0.05. 

Results 

Olfactory Test Results 

None of the subjects reported any perception of pathological 
pain or smell sensitivities. Most subjects (n = 70) were normosmic 
as indicated by a TDI ("Threshold Discrimination Identification") 

score >30.5 [14]. Neither the subjects' sex (x^ test: p = 0.45) nor 
their smoking habits were associated with the diagnosis of 
hyposmia (x^ test: p = 0.92). However, the TRPAl rsl 1988795 A 
allele was rarer among the 11 hyposmic subjects (TDI of 24.75- 
30.75). Specifically, the minor allele was carried by two (18.2%) of 
the 11 h^-posmic subjects and by 34 (48.6%) of the 70 normosmic 
subjects (cross-talnilation: likelihood ratio = 3.87, p = 0.049; 
Figure 1). 

Odor task results were diflferentiy affected by the TRPAl 
genotype (rm-ANOVA effect "odor task": df= 3,237, F = 0.041, 
p = 0.99, interaction "odor task" by "genotype": df= 2,237, 
F=3.34, p = 0.02; with a tendency toward significant between- 
subject effects of "genot>pe", df= 1,79, F = 3.213, p = 0.077). 
Additionally introducing "gender" into the analysis failed to 
produce significant effects (p>0.05 for all main effects and 
interactions involving "gender"). In the post hoc tests, better 
olfactory performance in carriers of the TRPAl rsl 1988795 A 
allele was observed with odor discrimination (t-test: p = 0.029; 
13.1±1.5 versus 12.3±1.6 correct discriminations, Cohen's 
d = 0.5). Of note, of the 16 odor discrimination items, butanol 
(target) could be discriminated from fenchone (non-target) by 35 of 
the 36 A allele carriers but by only 34 of the 45 wild type subjects, 
which was significant after a correction (cross-tabulation: likeli- 
hood ratio = 8.76, p = 0.003). The most pronounced genotype 
effect was observed on the perceived intensity of the H2S stimuli. 
Specifically, carriers of the A allele perceived the H2S stimuU 
significanfly more intense than non-carriers (29.2±13.2 versus 
21±12.8mm VAS, t-test: p = 0.006, Cohen's d = 0.63), which 
persisted after a correction (p = 0.024). The significance of this 
effect also persisted when excluding hyposmic subjects from the 
analysis (p = 0.044). Moreover, the ratings of the H2S stimulus 
intensity correlated with the perceived cold (Spearman's p = 0.28, 
p = 0.012, and p = 0.27). 
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Figure 1. Olfactory parameters and their modulations by the TRPA1 rs11 988795 G>A genotype (wild type subjects, GG, versus 
carriers of the minor A allele). From left to right: (i) In cross-tabulation of normosmic (N) versus hyposmic (H) subjects, hyposmic subjects were 
underrepresented among the carriers of the A allele (likelihood ratio test: p = 0.049). (ii) This genotypic association was more produced in the 
performance in the odor discrimination subtest of the Sniffn' Sticks test battery [18], where carriers of the A allele performed significantly better (t- 
test: p = 0.029). (ill and iv): The other olfactory tests, however, were unaffected by the TRPAl genotype (odor threshold and odor identification), (v) 
The genotype effects were most pronounced in the perceived intensity of olfactory H2S stimuli administered with an olfactometer (p = 0.006). The 
widths of the bars or boxes are proportional to the respective numbers of subjects per group. *: p<0.05, **: p<0.01. The minimum, quartiles, median 
(solid horizontal line within the box), and maximum are used to construct the box plots. 
doi:10.1371/journal.pone.0095592.g001 



Nociceptive Test Results 

In contrast to olfaction, nociception was not modulated by the 
TRPAl genotype in the whole cohort (p = 0.385-0.757 for the 
various nociceptive measures, Cohen's d = 0.07-0. 196). Therefore, 
genetic modulations were searched in subgroups of subjects. Four 
pain sensitivity subgroups were identified (Table 1 and Supporting 
Figure 1). The TRPAl variant exerted distinct genotype effects on 
different pain parameters (rm-ANOVA effect "pain measure": 
df= 5,365, F= 2.596, p = 0.025, interaction "pain measure" by 
"pain cluster": df= 15,365, F= 14.12, p<0.001, effect "geno- 
type": df= 1,73, F= 10.138, p = 0.001, interaction "genotype" by 
"pain cluster": df=3,73, F = 3.205, p = 0.018; further, non- 
significant, effects and interactions not shown). In analyses of 
single pain stimuli, cold (ANOVA factor "genotype": df= 1,73, 
F = 8.305, p = 0.005, Cohen's d = 0.53, 0.9, 0.12 and 2.5 in 
clusters #l-#4, respectively) and heat (df= 1,73, F = 9.272, 
p = 0.003; Cohen's = 1.4, 0.17, 0.48 and 3.2 in clusters #l-#4, 
respectively) pain thresholds were modulated by TRPAl 
rsl 1988795 OA (Figure 2). The modulation of heat sensitivity 
was not uniform across clusters (interaction "genotype" by "pain 
cluster": df=3,73, F= 7.056, p = 0.0003). Other pain readouts 
including the intensity of CO2 were unaffected by the TRPAl 
genotype (48. 2± 16. 8 mm VAS in A allele carriers versus 
51. 4± 15.6 mm VAS in non-carriers, t-test: p = 0.39). Additionally 
introducing "gender" into this analysis failed to produce 
significant effects (p>0.05 for all main effects and interactions 
involving "gender"). 



Discussion 

The TRPAl variant rsl 1988795A was associated with enhanced 
perception of odorous stimulants. The direction of this genotypic 
effect toward enhanced sensitivity agrees with the expectation 
from a previous observation, namely that the same variant 
enhanced sensitivity toward nociceptive cold stimuli [12]. This 
suggests that a functional association of the human TRPAl variant 
is not restricted to the perception of pain but extends to odour 
perception. However, other than with Na^,1.7 variantswhere 
effects of a gain-of-function variant in the sodium channel gene 
SCM9A [13] can be attributed to the olfactory system, with TRPAl 
a trigeminal component of the observed genetic modulation of 
odorant perception needs clearly to be considered. 

Several substances that evoke an olfactory sensation have been 
shown to excite TRPAl channels. This particularly applies to 
butanol for which the genotype effect on its discrimination from 
fenchone was presently observed. Butanol is a primary alcohol 
which at 1 mM applied to HEK293 cells expressing hTRPAl 
increased Ca^"*^ currents [32] . However, the present study was not 
sufficiently powered to observe patterns in the correct discrimina- 
tions of 16 target/non-target odor combinations and therefore, 
statistical findings have to be regarded with caution. A presence of 
TRPAl genotypic effects for TRPAl excitants cannot be claimed 
from the present observations. Several further odorants have been 
associated with TRPAl excitation such as eugenol [7] and 
amylacetate [11], however, this did not provide a clear pattern 
with respect to genotype differences in odor discrimination. In 
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Figure 2. Experiment to re-establish the previously observed association of the TRPA 7 rs11 988795 G>A variant with cold (and heat) 
pain [12]. A statistically significant genotype effect on heat and cold pain thresholds (ANOVA factor "genotype") was most pronounced in cluster #4 
characterized by a high heat pain sensitivity (Table 1), which was reflected in a significant ANOVA interaction "genotype" by "pain cluster" for heat 
pain thresholds. The minimum, quartiles, median (solid horizontal line within the box), and maximum are used to construct a "box and whisker" plot. 
doi:1 0.1 371 /journal.pone.0095592.g002 



addition to the insufficiency of the power, the relative excitation by 
butanol and other excitants presently cannot be judged since the 



observations originate from different studies [1 1,32]. As n-butanol 
had been used as test substance for olfactory thresholds in the first 
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Table 1. Decision rules (separated by lines) extracted from the CART classifier, providing a semantic description of the pain 
phenotypes (clusters #l-#4) found by Ward cluster analysis. 



Case belongs to 




IF (Rule conditions) 


Cluster #1: 


(IF 


CO2 intensity <42. 5 mm VAS 




AND) 


Heat threshold >46.55'C 




OR (IF 


42.5 mm VAS<C02 intensity <50.25 mm VAS 




AND) 


Pressure pain threshold >27.25 N/m^ 


Cluster #2: 


IF 


CO2 intensity >50.25 mm VAS 




AND 


Pressure pain threshold >27.25 N/m^ 


Cluster #3: 


IF 


CO2 intensity >42. 5 mm VAS 




AND 


Pressure pain threshold <27.25 N/m^ 


Cluster #4: 


IF 


CO2 intensity <42. 5 mm VAS 




AND 


Heat threshold <46.55=C 



The CART identifiers of the cluster membership correctly assigned 95.1% of the subjects. 
doi:l 0.1 371 /journal.pone.0095592.t001 



version of the Sniffii' Sticks test (e.g. [14]) before it had been 
replaced with phenyl ethyl ethanol, a possibility that the observed 
effects were related via the olfactory rather than the trigeminal 
system remains. However, while this would receive support by the 
expression of TRPAl in the olfactory epithelium and olfactory 
bulb in mice, being specific olfactory structures, the translation to 
humans remains uncertain. TRPAl is absent from the human 
olfactory bulb [33,34] and the expression pattern in the human 
olfactory epithelium is unknown. Therefore, concluding a purely 
olfactory effect in humans would stand on weak grounds. 
Moreover, the olfactory bulb is in close relation with trigeminal 
afTerents [35-37] which further impedes a clear distinction of the 
sensory qualities affected by the TRPAl variant. Finally, olfactory 
tests results are known to be influenced by cognitive factors 
[38,39], which had not been acquired presently. However, 
considering these factors as a potential confounder of the present 
results would imply that the TRPAl variant had an influence on 
them, for which no clear mechanistic hypothesis exists. 

The perception of the smell of terpineol, amylacetate, benzal- 
dehyde and toluene involves TRPAl excitation, however, they 
seem to be detected as trigeminal irritants [11]. Therefore, the 
presently observed changes in olfactory measures might have been 
conferred via a trigeminal mechanism. This would explain the lack 
of TRPAl' s genetic association with odor thresholds, which had 
been obtained using the purely olfactory stimulant phenylethy- 
lethanol. Apparently contradicting this hypothesis is the observed 
increase in the perception of H2S being a pure olfactory stimulus 
that cannot be perceived by anosmic subjects [40,41], thus 
providing that no unwanted trigeminal co-stimulation had taken 
place. To pursue the hypothesis of a trigeminal component, we 
performed an additional experiment. Specifically, seven subjects 
received 15 stimuli of clean air (200 ms, ISI 20-30 s) and 15 no- 
stimuli, i.e., no olfactometer action, at random order. The 
question, presented within 2-5 s after each stimulus, whether or 
not they had received a stimulus was correctly answered 30, 28, 
24, 24, 26, 17 and 18 times by subject 1-7, respectively, which was 
always more often than by pure chance (15 correct responses; 
test: p<0.0001). This indicated trigeminal co-stimulation and 
further strengthens the association of the observed effects with the 
trigeminal rather than olfactory sensory system. However, it was 
very slight and it was only detected when specifically addressed, 
not during the main experiments. A slight trigeminal co- 
stimulation is known to play a key role in the unconscious 



detection of odors [42] . As this effect has been described to occur 
more at a subconscious level [42], its absence on suprathreshold 
trigeminal CO2 stimuli might be explicable. Thus, the unintended 
trigeminal co-stimulation provides an explanation of the increased 
HgS intensity. It is better compatible with the absent effect on odor 
threshold than a direct TRPAl effect on olfaction, which however 
remains an alternative possibility when assuming that H2S is an 
excitant of TRPAl, in contrast to phenylethylethanol. This is 
conceivable as odorants differ in their effects on TRPAl channels 
[11]. However, a direct involvement of TRPAl in human 
olfaction is conceivable although current knowledge provides no 
clear proof. That is, TRPAl expression in the olfactory epithelium 
[8] and olfactory bulb [9] has been shown to occur in mice. In the 
human olfactory bulb, however, TRPAl was neither found by a 
proteomic analysis [33] nor by mRNA quantification on a 
microarray [43]. Whether it is expressed in the human primary 
olfactory neurons remains to be shown. 

The clustering analysis identified a subgroup of subjects sharing 
similar pain phenotypic patterns. On one of these patterns, the 
TRPAl variant could be associated with different sensitivity to 
thermal pain among carriers and non-carriers. However, the 
clustering resulted in a split of the sample into smaller subsamples 
which jeopardized the statistical power. On this basis, the present 
results regarding nociceptive signals have to be viewed with 
reservation although the genetic association with pain was 
implemented mainly to address compatibility with previous results 
[12]. It has to be noted that the reproduction of the mentioned 
results was not exact as the original association was observed in 
heterozygous subjects whereas in the present analysis, heterozy- 
gous and homozygous carriers of the variant TRPAl rsl 1988795 A 
were pooled to maintain statistical power. A nevertheless 
performed analysis with three genotype groups showed a 
significant between-subjects effect of "genotype" (df = 2,69, 
F = 8.26, p = 0.001). Moreover, contrasting heterozygous carriers 
with non-carriers while leaving out the homozygous carriers also 
showed a significant genotype effect (df= 1,66, F = 9.705, 
p = 0.003) stratified for pain clusters (interaction "genotype" by 
"pain cluster": df = 3,66, F = 3.853, p = 0.013). However, in the 
whole cohort no strong association with pain readouts could be 
observed. As we have recently demonstrated [44], clustering pain 
phenotypes provides larger genetic effect sizes, which in the 
present analysis had to be obtained with another commonly 
recommended [45] clustering strategy due to the small sample size. 
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Indeed, clustering strategies are increasingly pursued in pain 
research [29,45] as they fulfill the assumption that different pain 
phenotypes are based on different molecular pathomechanisms 
that are accessible to specific treatments. Reproduction of non- 
clustered data have often failed (e.g., [46] versus [47]), not 
necessarily because of non-functionality of the genetic variants. In 
a randomly chosen sample, the pattern of functional variants that 
are concomitandy present with the variant of interest can be also 
as such that a canceUing-out of the effects occurs [48] . This leads 
to small overall effects sizes [49] while greater effects apply to 
subgroups [50]. A consequents of tiu; need for clustering is that 
genetic association studies in pain require an increasing number of 
participants. The statistical power effect further decreases when 
gender becomes an additional factor modulating the genotype 
effect, which in a previous study, however, has been shown to not 
apply to the present TRPA] rsl 1988795 G>A variant [.30] but in 
the same study associations of variants in other genes were 
modulated by the subject's gender (e.g., variants in TRPVl, 
OPRDl, COMT). 

Interestingly, the present finding reproduced the previous 
observation of an association of the TRPAl rsl 1988795 variant 
with different sensitivities to thermal pain. Specifically, carriers of 
the minor A allele had a shorter withdrawal latency to noxious 
cold and in a small subgroup, the allele was also associated with 
differences in the sensitivity to heat pain. Indeed, role of TRPAl 
channels in thermal sensation is controversially discussed [51]. 
The observation of an association with heat pain had been 
previously interpreted on the basis of an earlier reported 
obser\'ation that local capsaicin administration affected the cold 
sensation longer than any other sensations including heat [52]. 
Those observations point at interactions among TRP channels 
[53]. TRPAl channels are often co-expressed with heat (>43°C 
[54]) gated TRPVl [4,55]) and the channels act in concert [56]. 
That is, TRPVl can oligermerize with other TRP family subunits 
including TRPV3 and TRPAl [57]. The heteromerization 
between TRPVl and TRPAl channels can substantially affect 
the calcium signaling pathways of TRPAl homomers [57] and the 
co-expression can lead to outward rectitication of single channel 
current-voltage relationships and modulation of open probabilities 
[58]. Moreover, while thermal (heat) hyperalgesia was initially 
attributed solely to TRPVl, currently TRPAl and TRPVl are 
regarded to be interdependently regulated downstream of PLC- 
coupled bradykinin (BK2) receptors [59] to establish hypersensi- 
tivity to heat [56]. 

It has finally to be noted that the molecular mechanism of 
functional associations of TRPAl rsl 1988795 is still unknown. The 
SNP is located in intron 20-21 of the TRPA] gene (Ensembl 
transcript ID ENST00000262209) with a distance of 618/908 
nucleotides from exon boundaries. In the HapMap [60] CEU 
cohort, it is located in a 25 kb haploblock spanning chromosome 8 
position 73102904 (rsl2550748) to 73128527 (rs3735942) (Sup- 
porting Figure 2). This haploblock contains two splice sites 
(rs3824151 and rs35427625) and one synonymous (rsl3280644 
G>A; L830L) variant. A TRPA] splice effect associated indirectly 
with rsl 1988795 may provide a hypothesis for the molecular 
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Conclusions 
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Supporting information 

Figure SI Heat plot and cluster dendrograms of the 
individual responses to pain stimuli. More intense orange 
coloring indicates higher pain sensitivity equal to lower pain 

thrc'sholds. Ward clustering suggested two or four pain sensitivity 
subgroups (Figure 2). As the latter singled out subjects with 
censored cold or cold/menthol thresholds (0°C) in a separate 
cluster, the four-cluster solution was preferred as it was not 
substantially worse that the two-cluster solution as assessed by 
means of Silhouette plots (not shown). 
(EPS) 

Figure S2 Haploblock structure and SNP localization of 
the TRPAl gene of the HapMap CEU population (Utah 
residents with Northern and Western European ances- 
try; http://hapmap.ncbi.nlm.nih.gov/). The TRPAl rsl 1988795 
G>A SNP (framed red) is located in a haploblock that comprises 
many more SNPs. Whili; rsl 1988795 is intronic: but not in a splice 
site, the haploblock contains 2 sphce site SNPs (rs3824151 and 
rs35427625) that according to this analysis are the most likely 
candidates to confer the repeatedly observed functional association 
of the SNP. Top: SPNs of the dbSNP database in the region of the 
identified TRPAl haploblock that contains rsl 1988795. Bottom: 
Haploblock structure according to a HaploView analysis, using the 
95% confidence intervals on D' [62] method as the default of the 
HaploView software [63] and showing the linkage between SNPs, 
those in the HapMap data set are given at the top of the 
haploblock structure. 
(EPS) 
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